Ultrafast Electron Dynamics Theory of Photo-excited Ruthenium Complexes 
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An explanation is provided for the ultrafast photo-excited electron dynamics in low-spin Ruthe- 
nium (11) organic complexes. The experimentally-observed singlet to triplet decay in the metal-to- 
ligand charge-transfer (MLCT) states contradicts the expectation that the system should oscillate 
between the singlet and triplet states in the presence of a large spin-orbit coupling and the absence 
of a significance change in metal-ligand bond length. This dilemma is solved with a novel quantum 
decay mechanism that causes a singlet to triplet decay in about 300 femtoseconds. The decay is 
mediated by the triplet metal-centered state ("^MC) state even though there is no direct coupling be- 
tween the ^MLCT and ^MC states. The coupling between the ^MLCT and ^MC via excited phonon 
states leads to vibrational cooling that allows the local system to dissipate the excess energy. In the 
relaxed state, the population of the "^MC state is low and the metal-ligand bond length is almost 
unchanged with respect to the initial photoexcited state, in agreement with experiment. 

PACS numbers; 78.47.J-, 33.50.-j, 82.50.-m 82.37. Vb 82.53.-k 
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Introduction.— Photon- induced nonequilibrium elec- 
tron dynamics in transition-metal complexes have drawn 
intense interest for numerous applications in diverse 
research fields such as solar energy conversion, pho- 
tosynthesis, photodissociation, and photocatalysis [ij. 
The study of the time evolution of the light-driven 
state is crucial in creating a deeper understanding of 
these technologically-important phenomena. Often, the 
photo-induced state undergoes an ultrafast decay into 
a metastable state from which it relaxes more slowly 
back into the ground state. A well-studied example is 
the spin-crossover phenomenon in Fe-complexes, where 
the spin-orbit coupling causes a cascade from a low-spin 
state to a high-spin state [2] ■ Crucial in the understand- 
ing of this process 0] is the change of the metal-ligand 
distance resulting from the transfer of electrons from 
the t2g to the Cg orbitals that repel the ligands more 
strongly jj]. However, a significant bond elongation is not 
present in all spin-crossover systems. In ruthenium com- 
plexes, such as [Ru''(byp)'^]^+, spin flips occur within the 
t2g orbitals and only a minimal change in metal-ligand 
bond length is observed Q. Due to their wide frequency 
response, ruthenium compounds have been extensively 
studied. The ground state is a singlet with six electrons 
in the t2g orbitals. Since Ru AdP multiplet excitations 
are dipole forbidden, the optical excitations are domi- 
nated by metal-to-ligand charge-transfer (MLCT) into 
the TT* orbitals of the organic ligands, creating a singlet 
^Igi configuration, where L indicates an electron in the 
TT* ligands. The strong 4d spin-orbit coupling then causes 
spin flips in the i^g states leading to triplet MLCT states. 
Since for the triplet t2gL state the coupling to the metal- 
centered (MC) and i^g^s configurations is either for- 
bidden or weak, it is lower in energy than the antibond- 
ing singlet states, see states 1 and 2 in Fig. [TJ The weak 
coupling has the interesting effect that it tends to local- 



ize the photo-excited electron on one of the bipyridine 
molecules [6]. The singlet-to-triplet decay is very fast. 
For example, in [Ru^^(byp)'^]^+, the relaxation time is 
300 femtoseconds (fs) [y, |7| . The resulting triplet state is 
long-lived with a lifetime on the order of hundreds of ns 
i|. 

Despite decades of research, the mechanism for the ul- 
trafast decay from the singlet to the triplet MLCT states 
is still not understood. The radiationless ultrafast decay 
is generally accompanied by a pronounced change of the 
metal-ligand bond lengths However, one does not 
expect a significant change in the Ru-ligand distance for 
different configurations. This has been confirmed ex- 
perimentally where for the ■^MLCT a bond contraction of 
only ~0.03 Angstrom (A) [5j] with respect to the ground 
state is observed. As a result, the Huang-Rhys factor 
(proportional to the square of the bond length change), 
and hence the decay, between MLCT singlet and triplet 
is expected to be negligibly small. Therefore, the sys- 
tem should oscillate between the singlet and triplet states 
via the spin-orbit interaction. Furthermore, it was al- 
ready noticed that a discussion in terms of singlet 
and triplet states is probably erroneous due to the large 
spin-orbit coupling with respect to the singlet-triplet en- 
ergy gap. However, this picture is completely contrary 
to the standard description of ultra-fast decay between 
spin-labeled states in Ru-complexes [1, 0] and, in almost 
all literature, a discussion of the mechanism underlying 
the singlet-to-triplet MLCT decay is generally avoided. 

In this Letter, we reconcile these apparent contradic- 
tions and present a coherent picture of the singlet-to- 
triplet decay in the presence of a strong Ad spin-orbit- 
coupling. Crucial in the mechanism is the presence of 
the triplet i^g^s metal-centered ('^MC) state, which has 
a Ru-N bond len gth about 0.4 A larger compared to the 
ground state ll|, . This state is generally considered 
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FIG. 1: Simplified energy level scheme of Ru^'''-based or- 
ganic complexes. The initial photon excitation from the sin- 
glet ground state (^GS) to the singlet metal-to-ligand charge 
transfer (^MLCT) state is followed by a relaxation to the long- 
lived metastable '^MLCT state. The triplet metal-centered 
^MC is located close to the MLCT states. 



to play an important role in the decay of the ^MLCT 
back to the ground state. There is consensus that the two 
triplet states are close to each other in energy l^, 13- 15|. 
The ^MC energy position with respect to the "^MLCT is 
experimentally controversial. However, the ■^MC is likely 
higher in energy than the ^MLCT, considering the sta- 
bility of the latter, see Fig. [TJ We demonstrate that 
the ^MC effectively mediates the decay from the ^MLCT 
to the '^MLCT even in the absence of a direct coupling 
between the ^MLCT and the ^MC state. The ^MC con- 
tribution to the relaxed state is minimal. 

Model and parameters.— In describing the singlet to 
triplet MLCT decay, we consider three excited states, 
iMLCT, 3MLCT, and ^MC, labeled as state 1, 2 and 
3, respectively. These levels couple to the vibra- 
tional/phonon modes of energy huj of the surrounding 
ligands, via the Hamiltonian 



Hq = EiUi + hua^ a + XiUi (a^ -I- a) 



(1) 



where rii gives the occupation of state i and is the 
creation operator for the vibrational mode. The cou- 
pling is dominated by the stretching mode that changes 
the metal-ligand bond length. Variations in coupling 
strength A,j lead to different metal-ligand equilibrium 
distances for the different states. The electron-phonon 
self-energy is Si = Xf/{huj). Since a translation of the 
coordinates shifts all the by a constant in Hq, only 
the relative change in coupling is of importance. There- 
fore, we define the electron-phonon self-energy difference 
Sij — (Ai — Aj)^ /{^^) between state i and j. In addition, 
the energy difference between the states after diagonal- 
ization of Hq is defined as A^- = {Ei — £i) — {Ej — Sj). 
The different states couple to each other with an inter- 
action of strength Vij , 



where cjc^ causes a particle-conserving transition be- 
tween states j and i. The size and nature of the coupling 
constants Vij is discussed below. The system is not sep- 
arated from its surrounding. Exchange of phonons with 
the bath damps the metal-ligand bond length oscillations 
via intramolecular energy redistribution processes. The 
vibrational relaxation by the bath is taken into account 
by the use of a dissipative Schrodinger equation [3, [3| • 
For a quantitative description, we need to determine 
the parameters in a reasonable range based on both ex- 
periment and theory. Fluorescence spectra show ^MLCT 
emission between 500 and 575 nm and '^MLCT features 
around 620 nm This corresponds to an energy dif- 
ference Ai2 = 0.3 eV. The electron-phonon self-energy 
difference e oc (Ax)^ is related to Ax, the difference of 
the metal-ligand equilibrium distances between two dif- 
ferent states. Since the singlet and triplet MLCT have 
almost the same metal-ligand equilibrium distances @, 
ei2 = 0.0 eV. The difference in metal-ligand bond length 
between the metal-centered and MLCT states is several 
tenths of an A, corresponding to a change of electron- 
phonon self-energy of the order of several tenths of an 
eV 0, Hi] and we take £13 = £23 = 0.4 eV. For the 
Ru-ligand stretching mode, we take a typical energy of 
huj = 30 meV ■ The coupling between the singlet and 
triplet MLCT states is the spin-orbit coupling of ruthe- 
nium, V12 = 0.12 eV [2^. There is no direct coupling 
between the ^MLCT and '^MC. In pure octahedral sym- 
metry, there is also no coupling between the '^MLCT and 
the ^MC, since the eg states would not couple to the tt* 
states. However, since most Ru-complexes have a lower 
symmetry, a small hybridization between the two triplet 
states should be present. We take a weak coupling con- 
stant of V23 = 0.08 eV, where the actual value will depend 
on the amount of distortion. The environmental relax- 
ation time for a single phonon mode is taken as 40 fs. 

Rate Equations.— First, we demonstrate that the de- 
cay from singlet to triplet MLCT states cannot be un- 
derstood from phenomenological rate equations based on 



the Franck-Condon factors [16|, |2l|, |22| . The relaxation 



is the Franck-Condon factor where n p 
g = 612/ huJ is the Huang- Rhys factor 



Hi = ^F,,(clc,-t-H.c.), 



(2) 



constant from state 1 (^MLCT) to state 2 (^MLCT) is 
r = 2TTFnV^2l^^^ where ^12 is the spin-orbit coupling 
between singlet and triplet levels and F„ = e~^.g"/7i! 

Ai2/?ia; and 
In the spin 

crossover in divalent Fe-complexes, electrons are trans- 
ferred from t2g to eg orbitals, resulting in a change of 
bond length on the order of tenths of an A. The corre- 
sponding e is tenths of an eV. However, in Ru complexes, 
the Ru-N bond length changes only about 0.03 A, i.e. an 
order of magnitude less Q . Correspondingly, the change 
in phonon self-energy £12 ~ (Aa;)^ should be two or- 
ders of magnitude smaller than in typical Fe-complexes. 
Since this is negligible with respect to A12, it is impos- 
sible to meet the condition for ultra-fast decay that the 
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FIG. 2: (color online) The probability of the ^MLCT (Pi(t), 
green thin solid), the ^MLCT (P2(t), red solid), and the ^MC 
(P3(t), blue dashed) as a function of time for different values 
of the energy difference A32 between the ^MC and ^MLCT 
states, (a) A32 = 0.35 eV; (b) A32 = 0.05 eV; (c) A32 = 
-0.15 eV. 



energy gap is near the electron-phonon selfenergy differ- 
ence The hfetime is therefore expected to be much 
longer than picoseconds (ps) . Since there is also no direct 
coupling between state 1 and 3 ('^MC), the initial state 1 
cannot reach state 2 either directly or via state 3. There- 
fore, phenomenological rate equations cannot provide a 
satisfactory explanation of the singlet to triplet MLCT 
ultrafast decay. 

Singlet-to-triplet MLCT decay— In the following we 
investigate the time-dependence of the states based on 
quantum- mechanical principles First, in the ab- 

sence of state 3 ("^MC), we find, as expected, no decay 
but an oscillatory behavior involving the triplet and sin- 
glet MLCT states as a result of the spin-orbit coupling 
(not shown). When the ^MC is above state 1 (^MLCT), 
the situation is very similar and no decay is observed, 
see Fig. da) for A32 = 0.35 eV (A31 = 0.05). The 
population of the ■^MC state is almost negligible. There 
is no decay and an an oscillating population of the sin- 
glet and triplet MLCT states. The oscillation period of 
11 fs is close to frequency expected from the energy gap 
a/ + 41^2 in tii6 absence of the ^MC state. 

Let us now consider the situation where the '^MC state 
is located between the MLCT states, see Fig. Eljb) with 
A32 = 0.05 eV (A31 = -0.25 eV). Ultrafast decay oc- 
curs, and after about 300 fs the average populations 
of the system become stable. The population of the 
^MLCT state increases to about 85% at 600 fs. Al- 
though, at first sight, the time dependence of the prob- 
abilities resembles a cascading from state 1 (^MLCT) to 
2 (^MLCT) via state 3 (^MC), it is important to remem- 
ber that in the model there is no direct coupling between 
state 1 and 3 (V13 = 0). Therefore, in a classical rate 
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FIG. 3: (color online) The metal-ligand bond length as a func- 
tion of time t for different energy differences A32 between the 
''MC and ^MLCT states. The bond length change is normal- 
ized to the bond elongation (xg = 0.4 A) of the MC triplet 
state with respect to the MLCT state The oscillation period 
is close to the phonon's Tphonon =137.8 fs. 

equation, the expected decay time would be very large, 
Tis = (2'KFnVi^/h'^u!)~^ — > 00. A detailed analysis of the 
evolution of the state occupations with time, shows that 
the initial ^MLCT state couples to the lowest vibrational 
level of the '^MLCT via the spin-orbit coupling. Minor 
distortions from octahedral symmetry allow a hopping 
V23 of the electron from the ligand tt* orbitals into the 
ruthenium Cg orbital, populating the '^MC state. How- 
ever, due to the large change in electron-phonon self- 
energy £23, the hopping to the lowest vibrational '^MC 
state is strongly reduced due to the Franck-Condon fac- 
tors. Therefore, the hopping between the ^MLCT and 
■^MC does not conserve phonons. This allows the sys- 
tem to decrease its energy via vibrational cooling due 
to intramolecular energy redistributions The prob- 
ability in the ■^MC state can however also couple back 
again to the ■^MLCT, again leading to the creation of ex- 
cited phonon states that are subject to vibrational cool- 
ing. Therefore, the interplay between the ■^MLCT and 
the "^MC state leading to the creation of excited phonons 
allows the local system to dissipate energy via the damp- 
ing of the Ru-ligand bond-length oscillations. The con- 
sequence is a reduction of the recurrence of the ^MLCT 
state. In the relaxed state (> 500 fs), the contributions 
of the '^MC is small after reaching a maximum of close to 
50% around 100-150 fs. The contribution of the ^MLCT 
is of the order of 10%, justifying the traditional designa- 
tion of the final state of the decay in Ru-complexes as a 
triplet. 

In addition, let us consider the case that the '^MC 
state is lower in energy than the ^MLCT state. For 
A32 = —0.15 eV, we observe a cascading process from 
the iMLCT via the ^MLCT to the ^MC state, see Fig. 
Wic). The '^MC reaches an occupation of approximately 
92% at about 600 fs. In a real system, the occupation 
of the '^MC state most likely decays back to the ground 
state (not included in the calculation). 

Finally, we study the evolution of the metal-ligand 



4 



bond length as a function of time, see Fig. [3l EXAFS 
measurements taken 70 ps after the photo-excitation 
show a bond length contraction of about 0.03 A with 
respect to measurements in the absence of a photoex- 
citation. This rules out a strong contribution of the 
■^MC, which is expected to have a significantly larger 
bond length due to the presence of electron density in 
the Cg orbitals. The small change in bond length is sup- 
ported by our model. The calculated bond length elon- 
gation is strongly related to the probability of state 3, 
see Fig. [5] For the situation, where the ^MC state is 
located between the MLCT state (A32 = 0.05 eV), the 
bond length initially increases and reaches a maximum of 
about 0.16 A around 100 fs. Subsequently, the length de- 
creases as the occupation of the ^MLCT state increases. 
At about 600 fs, the average bond elongation is only 
0.024 A larger with respect to the MLCT states. This 
should be contrasted with the situation where the ^MC 
state is lower in energy than the ^MLCT (A32 = —0.15 
eV). Here, the bond length increases by several tenths 
of an A, as the ^MC is occupied. On the other hand, if 
the ^MC state energy is higher than that of the ^MLCT 
(A32 = 0.35 eV), the bond length is also almost the same 
as that of the MLCT states. 

Conclusion.— Experiments on ruthenates over several 
decades have established that the photo-excited singlet 
metal-to-ligand state decays into its triplet counterpart 
with very little change in lattice parameter |^, |23| . 
This contradicts the expectation that, in the absence of 
a clear decay mechanism, an oscillation should occur be- 
tween the singlet and triplet MLCT states as a result 
of the strong ruthenium spin-orbit coupling. This paper 
has shown that a triplet metal-centered state located be- 
tween the MLCT states can effectively mediate the decay 
from the MLCT singlet to the triplet state. The decay 
cannot be simply described in terms of a phenomenolog- 
ical rate equation, but is a result of vibrational cooling 
of the oscillations of the Ru-ligand bond length due to 
the intersystem crossings between the ^MLCT and the 
•^MC states. The presence of the '^MC also strongly re- 
duces the mixing of singlet and triplet states justifying 
the standard singlet and triplet notation which at first 
appears erroneous for a compound with a large spin-orbit 
coupling. The decay occurs on a timescale of 300-400 fs 
and has a very high quantum efficiency. 

Several aspects provide interesting experimental tests 
of the described mechanism. First, although the 
metastable state reached after several hundreds of fem- 
toseconds shows little change in the metal-ligand bond 
length, there is an increase in bond length of the order of 
0.15-0.2 A in the 200-300 fs following the photoexcitation. 
The observation of this dynamic response is in the reach 
of state-of-the-art optical pump/x-ray probe techniques 
[3|. Second, the coupling between the ^MLCT and ^MC 
states is a result of finite hopping matrix elements due 
to distortions that break the octahedral symmetry. The 



dependence of the decay time in different geometries can 
provide additional information on the detailed nature of 
the decay mechanism. 
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